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Effective action of center vortices in SU (2) lattice gauge theory is investigated by studying the 
correlation between the action density on their worldsheets and their geometric properties. It 
turns out that center vortices are rigid, however, their dynamics is more complicated than that of 
rigid random surfaces, since some coupling constants have nonstandard scaling dimensions. As a 
result, the properties of center vortices are almost completely determined by curvature-dependent 
effects. This, in turn, provides a qualitative explanation of vortex percolation. 



The XXV International Symposium on Lattice Field Theory 
July 30 - August 4 2007 
Regensburg, Germany 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Rigidity and percolation of center vortices 



M. I. Polikarpov 



1. Introduction 

Since the seminal work of t'Hooft [1] center vortices in non-Abelian gauge theories are often 
discussed as field configurations which are responsible for confinement. During the last decade the 
properties of center vortices were extensively studied using lattice simulations, and their impor- 
tance for the low-energy behavior of Yang-Mills theories has been fully verified: it was found that 
confinement and chiral symmetry breaking are absent if vacuum expectation values are calculated 
for lattice configurations without center vortices [2]. In [2, 3] it was demonstrated that the area of 
center vortices in SU (2) lattice gauge theory scales in physical units of length, which gives one 
more argument that center vortices are not just lattice artifacts. Many other interesting proper- 
ties of center vortices were discovered in [3, 4]. For instance, it was found that center vortices are 
genuinely thin and carry quadratically divergent action density. It also turned out that vortex world- 
sheets are densely covered with AbeUan monopoles [4]. However, a proper theoretical description 
of vortex dynamics is still to be found. 

In recent works [3, 5] an attempt has been made to find an effective vortex action by studying 
the correlation between the action density on vortex worldsheets and their local geometry. Center 
vortices turned out to be rigid, i.e. extra energy is required to bend the vortex. It was conjectured 
that rigidity may be induced by some two-dimensional fields locaUzed on vortices. Rigidity of cen- 
ter vortices allows one to understand at least quaUtatively the observed scaling of vortex area, since 
in contrast to the simplest case of random surfaces with Nambu-Goto action rigid random surfaces 
may have a nontrivial continuum limit [6]. In fact rigidity is a common feature of any vortex-like 
solution in gauge theories, since the terms in the string action which make vortex worldsheets rigid 
are the first terms in the general gradient expansion of the effective vortex action [7]. 

On the other hand, there are some important properties of center vortices which are not so 
easily explained by the model of rigid strings. One of them is the percolation of vortex world- 
sheets, which can play a crucial role in the confinement mechanism (for instance, it is known that 
percolating vortex disappears at the deconfinement phase transition [8]). The aim of this paper is 
to investigate the effective vortex action in more details than in [5] and to discuss how such action 
can describe the emergence of percolating vortex clusters. 

2. Effective vortex action 

A basic assumption of [3, 5] and this work is that the average total excess of action on the 
vortex worldsheets is equal to the effective vortex action Wgff [E] obtained by integrating over all 
degrees of freedom in the theory except the geometry of the vortex worldsheet E. In general, 
We//[r] can include nonlocal terms, for instance, of the form R^^R, where R is the internal 
curvature of the vortex worldsheet and A is the covariant two-dimensional Laplacian. As it is very 
difficult to trace such nonlocal terms on presently available lattices, in this paper we will assume 
that Weff [E] is indeed local. Furthermore, we will simplify our task by assuming that Wg// [E] 
depends only on the simplest local two-dimensional geometric invariants - internal curvature R and 

extrinsic curvature K, i.e. Wgff [S] = / (f'^ sf^^eff are the coordinates 

z 

on the vortex worldsheet, gab = daX^dbX^ is the induced metric and g is its determinant. The 
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extrinsic curvature K can be written in the following form: 

K = g'^^dj^'^dbt^" = -t^'^M^" (2. 1) 

where t^^ = e"'^ daX^dhX^ is the two-form tangent to the vortex worldsheet Z. We has also taken into 
account that the effective action should depend on the lattice spacing a, according to the Wilson's 
renormalization group equations. 




Figure 1: Average excess of action as the function of extrinsic (left) and internal (right) curvatures and 
lattice spacing 

In order to study the correlation between the excess of action density and internal and extrin- 
sic curvatures of vortex worldsheets in pure SU (2) Yang-Mills theory, the excess of action per 
plaquette ^Sp = a^We/f {K,R;a) = j8 (l - ^TrUp) - (jS (l - ^TrUp)) = f {{TrUp) -TrUp) was 
averaged separately for vortex plaquettes and sites with different values of K and R. Lattice dis- 
cretization of ( |2.1[ ) can be written as: 

a% = -t^^a'At^p- = l^t^^ (t^^ - =4 -Up (2.2) 

p' 

where ^ denotes summation over all plaquettes p' adjacent to the plaquette p, rip is the number of 

p' 

vortex plaquettes which are adjacent to p and are parallel to it and the factor comes from the 
definition of derivatives on the lattice. According to (^^, a^Kp takes integer values from to 4. 
Internal curvature for lattice surfaces was defined as in [5, 6] : a^R^ = 4 — n^, where is the number 
of sites on the vortex worldsheet which are adjacent to s. Thus in contrast to [5], where both K and 
R were associated with lattice sites, here K is associated with vortex plaquettes and R - with lattice 
sites on vortex worldsheets, which yields more self-consistent results. 

Average excess of action per vortex plaquette ( ASp ) as the function of lattice spacing and 
extrinsic and internal curvatures (in lattice units) is plotted on Fig. |l|. Empirically it was found that 
Weff {K,R;a) is almost linear in o'^Kp and has a distinct minimum with respect to o'^Rg, which is 
situated near a^R^ ^ —I. For the sake of simplicity we have approximated the dependence of ASp 
on a^Rs by a quadratic function, since in fact only the existence of this minimum is essential. The 
actual dependence of A^^^ on a^R^ at different lattice spacings and the quadratic functions which 
approximate it are plotted on Fig. ^ at lower right. Thus we assume that Weff {K,R;a) has the 
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following form: 

Weff{K,R;a) = o{a) + K{a)K + 'Y{a) {R-Ro{a)f (2.3) 

where Ro{a) ^ — I- a^^. This action differs from the action of rigid strings used in [5] by the terms 
quadratic in R. Although the standard dimensional analysis shows that only a {a) and the terms 
Unear in K and R should survive in the continuum limit, we shall see that for center vortices higher- 
order terms can be equally important. It is also more convenient to redefine a {a) = oT^Oq (a) = 
A^yOQ (a). The coefficients Oq (a), K{a) and /(a) as the functions of lattice spacing are plotted on 
Fig. I 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 -5 -4 -3 -2 -1 1 2 

a, fm Rp 

Figure 2: At upper left: bare vortex tension Go (a), at upper right: coupling constant K{a), at lower left: 
coupling constant 7(0), at lower right: average excess of action per plaquette as the function of internal 
curvature in lattice units a^Rp at different values of lattice spacing (from a = 0.14/m for the uppermost 
curve to a = 0.06 fm for the lowest curve. The errors in experimental results are not shown in order not to 
clutter the plot). 



3. Scaling of coupling constants and renormalization of the effective vortex action 

Dimensional analysis shows that as the continuum limit t? ^ is approached, the bare vortex 
tension a (a) should grow as a^^ ~ A^y and the coupling constant 7(a) should tend to zero as 
a^. Since the external curvature is a marginal operator, the coupling constant K"(a) can grow not 
faster than In a. 
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Quadratic divergence of surface action is a necessary condition for the physical scaling of the 
area of random surfaces. The reason is that for self-avoiding surfaces the number of surfaces of a 
given area S grows as exp (c A^yS) , where c is some constant which depends on the structure of 
the lattice, and the sum over all surfaces can only converge if the statistical weight of each surface 
decreases not slower than exp (— cA^^^). In [3] it was shown that the average action associated 
with center vortices is indeed proportional to the total area of the vortex worldsheets and diverges 
quadratically: 

{Weffil.)) ^OMa-^ S[Z]=OM^uyS[I.] (3.1) 

Lattice simulations indicate that Oo{a) is almost hnear in a (see Fig. ^: ao(a) ~ A + Ba, where 
A = 0.05 ±0.02 and B = (2.9 ±0.02) fm-\ therefore a (a) = Aa'^ +Ba-^. The value of A is 
numerically very small, and it is Ukely that a {a) diverges only linearly. Thus the bare vortex 
tension a {a) diverges too weakly to saturate ( |3.l| ) and to cancel the quadratically divergent surface 
entropy, hence quadratic divergence (|3.1|) in the effective vortex action should arise due to some 



other terms. If this action indeed depends on K and R only, the fit ( |2.3D and the estimation ( p.lD 
imply that curvature-dependent terms should diverge, i.e. K{a) {K) + y{a) { {R — Rq (a))^ ) ~ A^y. 
In order to check this we have measured the average number of parallel adjacent plaquettes {np ) 
and the average number of neighbors per site ( « v ) for vortex worldsheets as the functions of lattice 
spacing (see Fig. Extrapolation to the continuum limit shows that both a^Kp = 4 — np and 

— a^Rs = ns — 4 are finite and positive as a — > 0, therefore ( K ) and ( R ) indeed diverge as and 

— a~^. As we will see further, it is important that ( /? ) is negative. 

Consider now the coupling constants K{a) and y{a). Extrapolation of the results plotted on 
Fig. ^to a = shows that K{a) is finite in the continuum limit: K"(0) = 0.18 ±0.02, while y{a) 
is proportional to a: y{a) = Ca = CAyy, where C = (0.0022 ± 0.0001) /m. Standard dimen- 
sional analysis shows that y{a) should be quadratic in a, thus both a {a) and y{a) seem to have 
nonstandard scaling dimensions: 1 and —1 instead of 2 and —2. 

Extrapolation of {a^Kp) =4 — {np) shows that ( A') « (3.45 ±0.01) a^^ as a ^ 0. According 



to (23), the contribution of (A') to the average action density is ( Swgff {K,R;a) ) = K{a) (K) ^ 
(0.18 ±0.02) (3.45 ±0.01) a^^ ^ (0.6 ±0.1) a^^, which agrees well with the estimation (|1|). It 
appears that the dominating contribution to the effective vortex action We/f [I] comes from the term 
K{a)K, therefore center vortices can be effectively described as rigid tensionless strings. 

As y{a) has scaling dimension —1, quadratically divergent internal curvature yields a cubic 
divergence in Weff{K,R;a) unless (/?) =R(){a). Extrapolation of {a^Rs) to the continuum limit 
shows that (/? ) is indeed close to Rq (a) if the averaging is performed only over percolating vortex 
(see Fig. |3|). However, when the averaging is performed over all vortices, (R) — Ro{a) again 
diverges quadratically, thus formally giving rise to a cubic divergence in We/f [£] . This fact implies 
that either our measurements are flawed by some lattice artifacts (such as incomplete separation of 
percolating and satellite vortices) or that even higher powers of R survive in the continuum limit 
and, consequently, our effective vortex action is not renormalizable. 

4. Percolation of center vortices 

By virtue of the Gauss-Bonnet identity ( / d^^ ^R ) ~ ( S [£])(/?) ~ 2 - 2( ^ [£] ), where 
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Figure 3: Left: average number of adjacent parallel plaquettes as the function of lattice spacing. Right: 
average number of neighbors per site as the function of lattice spacing. 



g [E] is the genus of the surface Z. Therefore if {R) is large and negative, the genus of vortex 
worldsheets is also large, i.e. there are a lot of handles on it. Since center vortices have finite 
density in the continuum limit, ( 5 [£] ) and the genus are proportional to the total volume 

of physical space. This agrees well with the results of lattice simulations, which indicate that the 
genus of the percolating vortex is indeed very large and grows linearly with the total volume of the 
lattice. Thus the action of the form ( |2.3[ ) provides also a rough qualitative explanation of vortex 
percolation: due to the term 7(a) (R — RQ{a))^, where Rq < 0, it is more advantageous to have 
a single connected surface with a large number of handles, so that {R) Rq < 0, than a large 
number of disconnected surfaces with (R) > 0. It is interesting to note that the effective vortex 
action of the form Weff {K,R;a) = a (a) + K{a) K + "/ {a)R with / (a) > [5] also favors negative 
R, which, however, cannot take arbitrarily large absolute values because of the positively defined 
term K{a)K. The decreasing of (/?) should stop at some negative value ~/?o('^)> which is UV 
divergent, because both K (a) and / (a) are dimensionless and therefore neither K nor R know about 



the A.QCD scale. Presumably the term y{a) {R — Rq (a)) in (23) is just an effective description of 
this phenomena. 

Since (K) and (R) diverge as Ajjy = a~^, a typical size of geometric structures associated with 
center vortices should be of order of several lattice spacings. Indeed, a typical vortex configuration 
on the lattice, besides a single percolating vortex, contains also a large number of small satellite 
vortices, which consist of only several lattice plaquettes. However, such small vortices can give 
only positive contribution to (/?), and negative UV divergent values of {R) can only be associated 
with percolating vortex. This implies that a typical length scale associated with percolating vortex 
should also be comparable with lattice spacing, which suggests that percolating vortex emerges 
as a result of purely statistical coalescence of small satellite vortices, just as in the case of usual 
percolation. On the other hand, it is difficult to explain the physical scaling of total vortex area by 
such purely statistical mechanism which works at distances of order of lattice spacing. It seems 
that tensionless rigid random surfaces and percolation theory are two complementary languages 
which describe one and the same phenomenon. 
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5. Conclusions 

The analysis performed above shows that the dynamics of center vortices is more complicated 
than that of rigid random surfaces. In particular, the bare vortex tension diverges too weakly to 
compensate the divergence of surface entropy, and the fluctuations of surface geometry can only 
be suppressed by the terms which depend on the curvature of vortex worldsheets. As a result, 
center vortices form some irregular geometric structure with characteristic size of order of several 
lattice spacings, so that {K) ~ A^y and (i?) ~ — A^y. This structure, however, is completely 
different from branched polymers which are a common problem for random surfaces with Nambu- 
Goto action [6]. In fact negative {R) indicates that the genus of vortex worldsheets grows with 
their area, which is a necessary condition for percolation. In contrast, branched polymers are 
usually assumed to have fixed topology [6]. It is interesting to note that two commonly assumed 
conditions of surface percolation, namely, vanishing of vortex tension and the existence of "false" 
vacuum [8] are fulfilled in our case, although in a rather odd way. Percolating vortex in some sense 
corresponds to the "true" vacuum with nonvanishing negative {R), while small satellite vortices 
can be considered as perturbative excitations over this vacuum. 

This work was partially supported by grants RFBR-DFG 06-02-04010 and DFG-RFBR 436 
RUS 1 13/739/2, RFBR 05-02-16306, and RFBR 04-02-16079, and by the EU Integrated Infras- 
tructure Initiative Hadron Physics (I3HP) under contract Rn3-CT-2004-506078. 

References 

[1] G. t' Hooft, On the phase transition towards permanent quark confinement. Nuclear Physics 
B 138 (1978) 1 - 25. 

[2] L. Del Debbio, M. Faber, J. Giedt, J. Greensite, and S. Olejnik, Detection of center vortices in 
the lattice Yang-Mills vacuum. Physical Review D 58 (1998) 094501. 

[3] F. V. Gubarev, A. V. Kovalenko, M. I. Polikarpov, S. N. Syritsyn, and V. I. Zakharov, Fine 
tuned vortices in lattice SU(2) gluodynamics. Physics Letters B 574 (2003) 136 - 140. 

[4] V. G. Bornyakov, P. Y. Boyko, M. I. Polikarpov, and V. I. Zakharov, Monopole clusters at 
short and large distances. Nuclear Physics B 672 (2003) 222 - 238. 

[5] P. V. Buividovich and M. I. Pohkarpov, "Center vortices as rigid strings." ArXiv:0705.3745 
[hep-lat], 2007. 

[6] J. Ambj0m, "Quantization of geometry." Lectures presented at the 1994 Les Houches 
Summer School, 1994. 

[7] M. Engelhardt and H. Reinhardt, Center projection vortices in continuum Yang-Mills theory. 
Nuclear Physics B 567 (2000) 249. 

[8] M. N. Chemodub and V. I. Zakharov, Magnetic component of Yang-Mills plasma. Physical 
Review Letters 98 (2007) 082002. 



7 



